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Abstract 

We determine some points on the finite size scaling curve for the 
correlation length in the two dimensional 0(3) and icosahedron spin 
models. The Monte Carlo data are consistent with the two models 
possessing the same continuum limit. The data also suggest that the 
continuum scaling curve lies above the estimate of Kim Q and Carac- 
ciolo et al [Q] and thus leads to larger thermodynamic values of of the 
correlation length than previously reported. 



In 1993 Kim JT| proposed using finite size scaling to obtain the thermody- 
namic value of the correlation length £ in the two dimensional {2D) nonlinear 
a model at given inverse temperature j3. He claimed that he could predict 
correctly (less than 2% error) values as large as 15000 from measurements 
taken on lattices not larger than 283 x 283. Kim's method was adopted by 
Caracciolo et al 0, who, after refining it and collecting data at a 180 pairs of 
L and j3 values claimed to have such control over the scaling curve that they 
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could predict values of the correlation length as large as 10 5 even though 
their largest lattice was only 512 x 512. 

Of course at fixed x = £(L)/L the value of £(2L)/£(L) depends upon L 
and the continuum value is the limit of this ratio as L — > oo. Both Kim and 
Caracciolo et al claimed that within their statistical accuracy, they could not 
detect any systematic drift of £(2L)/£(L) with L within the range of values 
of L they studied; so they concluded that they had reached the continuum 
limit. (The latter authors also investigated in more detail the approach to the 
continuum limit in a different model, the two-dimensional SU(3) principal 
chiral model || and convinced themselves that their continuum extrapolation 
was supported by their data (private communication from Alan Sokal)). 

Shortly after the appearance of these papers, we criticized both Kim's 
paper [|J and the Caracciolo et al paper || by pointing out that we know 
rigorously that in the true continuum limit large spin fluctuations have 
to occur for any value of x > 0. Indeed it can be shown that the spins must 
become decorrelated over distances that are fixed, nonvanishing multiples 
of the infinite volume correlation length, while in the regime in which these 
authors took their data, at least at x > 0.75, the lattices were so well ordered 
that second order lattice perturbation theory (PT) could reproduce their 
Monte Carlo results within the errors. 

In this paper we significantly improve Caracciolo et al's statistics, and go 
to larger lattices (their largest L-2L pair has L = 256, ours 640). Thereby we 
find (for £(L)/L = 0.5) some statistically significant corrections to scaling 
even at the largest L values, which in our view signifies that the data of 
Caracciolo et al were rather far from the continuum limit even at x — 0.5 (of 
course, it is to be expected, and the present study corroborates that, that 
the continuum limit for larger values of x is reached at larger values of L). 

Our interest in this problem was revived by the recent discussions of 
symmetry enhancement of discrete nonabelian groups. Namely both we |7j || 
and Hasenfratz and Niedermayer |J pointed out that numerics suggest that 
the dodecahedron and icosahedron spin models have the same continuum 
limit as the 0(3) model. Hasenfratz and Niedermayer interpreted this as 
meaning that these discrete spin models were asymptoticaly free (AF). In fact 
in our paper we P| had measured the Liischer-Weisz- Wolff running coupling 
constant |10| and showed that it did not vanish at the critical point, as 



required by AF. A theoretical argument was also published shortly thereafter 
fllTf claiming to show the impossibility of the Hasefratz-Niedermayer scenario 
of AF in a discrete spin model. 
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This latter paper contains an interesting suggestion: the observed agree- 
ment between the discrete spin models and 0(3) is a transient phenomenon, 
which should disappear at larger correlation length. The observation stems 
from the fact that if one accepts the standard scenario of 0(3) being AF, then 
the discrete spin models could be regarded as a perturbed 0(3) model, and 
in the accepted PT scheme around the Gaussian fixed point, this would be 
a relevant perturbation, hence the two models could not possibly be equiva- 
lent. However the authors find that this perturbation becomes relevant only 
for (3 sufficiently large. They translate their estimate of this value of (3 into 
a correlation length of approximately 200 and suggest that perhaps the nu- 
merics will show that as the correlation length is increased, for continuum 
observables, the agreement between 0(3) and the icosahedron or dodecahe- 
dron spin models improves for a while, then, as £ exceeds approximately 200, 
it starts deteriorating. 

We decided to investigate numerically this possibility by comparing the 
finite size scaling curve of £(L) in the two models. We studied the icosahedron 
spin model and the 0(3) model at x — 0.25, x = 0.5 and x = 0.75. For the 
icosahedron we also took data at our estimated (3 crt = 1.803 ||. The results 
are recorded in Tabs. 1-7 and displayed in the figures. 

Before discussing the figures, let us specify in detail what we did and sum- 
marize what we find. Our systems consist of two dimensional square arrays of 
spins of size L x L with periodic boundary conditions. The spin at each site 
is of unit length and takes values either on the sphere S 2 or on the vertices of 
an inscribed regular icosahedron. Each spin interacts ferromagnetically only 
with its 4 nearest neighbours at inverse temperature f3. 

We used the same definition of the correlation length £ as 0: let P = 
(p, 0), p = n = 0, 1, 2, L — l. Then 

where 

G(p) = ^(|S(P)| 2 ); S(P) = £e^) (2) 

At a given L, we adjusted (3 so that the ratio £(L)/L was approximately 
0.25, 0.5 respectively 0.75 (we allowed differences only smaller than 2 x 10~ 3 ). 
Leaving (3 unchanged, we then doubled L and measured £(2L) and therefore 
£(2L)/£(L) at our x. (We used the slope of the scaling curve of || to correct 
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for the fact that our £(L)/L is not exactly equal to the desired values of 
0.25, 0,.5 and 0.75 ; since this correction is tiny, it does not matter whether 
their scaling curve represents exactly the truth). In our study L was varied 
from 20 to 640 or to 1280 and the data were produced using the usual one 
cluster algorithm. For each (3 L pair we performed several runs (up to 200), 
one run consisting of 100,000 clusters used for thermalization and 1,000,000 
for taking measurements. Each run started from a new randomly chosen 
configuration. The errors were computed from the results produced by the 
different runs, using the jack-knife method. 

In Fig.|I]we show our results for 0(3) (full squares) and for the icosahedron 
(open circles) at x — 0.5. For comparison, we also plotted the results of 
Caracciolo et al for 0(3) (full triangles). The choice of the abscissa 1/ (lnL + 
c) , was motivated by the fact that in studying lattice artefacts for the LWW 
step scaling function we found that such an ansatz seemed to describe the 
data pretty well (see ref.[|L2| Fig. 7). The value of the parameter c = 0.7 was 
obtained from a joint fit to our 0(3) and icosahedron data with a common 
limiting value for L — > oo; it is of acceptable quality (x 2 /dof = 11.6/8). The 
solid curve is a Symanzik type fit (a + b/L 2 + c\iaL/L 2 ) to the 0(3) data 
with, however, an unacceptable x 2 /dof = 10.7/3. Several facts are suggested 
by this figure: 



• Our MC data agree very well with those of Caracciolo et al. However 
while the largest L value investigated by these authors was L = 256, 
ours is L = 640. Also our error bars are much smaller. 

• The 0(3) data suggest a systematic increase of £(2L)/£(L) with L. This 
is true about both our data and those of Caracciolo et al, the latter 
however, having much larger error bars, could also be interpreted as 
showing no L dependance for L > 64. 

• The data suggest that the continuum value is definitely larger than the 
value predicted by Caracciolo et al 1.5420(7). 

• There is no obvious reason to suspect that the continuum limit in the 
two models is different. 

• The dashed lines represent linear extrapolations in l/(lnL + 0.7) and 
they intersect at £(2L)/£(L) = 1.584. The extrapolations used and the 
value they produce should be regarded only as an illustration, since we 
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Figure 1: Scaling function £(2L)/£(L) versus l/(lnL + 0.7) for 0(3) (full 
triangles: ref. and full squares : our data) and icosahedron (open circles) 
at £(L)/L = 0.5. The straight lines are fits. 



have no way of knowing whether our L values are truly asymptotic and 
the extrapolation ansatz is correct. 

For a better test of the approach to the continuum, we investigated a 
smaller value of x, namely x = 0.25, where one would expect the continuum 
limit to be reached at smaller L values. The data are shown in Fig. |2|. 

They suggest the following facts: 

• Indeed the continuum limit seems to be reached at lower L values, 
perhaps as low as L = 320 corresponding to a correlation length £ = 80. 

• The data are consistent with the two models sharing a common contin- 
uum limit value for the step scaling function, item This value seems 
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Figure 2: Scaling function £(2L)/£(L) versus l/(lnL + 0.7) for 0(3) (full 
squares) and icosahedron (open circles) at = 0.25 



to be around 1.031, higher than the value quoted by Caracciolo et al 
(1.0255(20)). 

In Fig.^ we present the same results as in in Fig.|l| but at x — 0.75. To 
go to this increased value of x, we had to increase j3. Consequently it is to 
be expected that at this x value, asymptopia will set in at a much larger 
value of L. This must be so especially for the 0(3) model, which will be 
described by PT up to much larger values of L (as we emphsized above, in 
the true continuum limit, the system cannot be in a PT regime, since in the 
continuum limit the spins decorrelate over distances proportional to L) . Two 
things can nevertheless be learned from this figure: 

• There is no reason to rule out that the continuum value is the same in 
the two models. 
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Figure 3: Scaling function £(2L)/£(L) versus l/(lnL + 0.7) for 0(3) (full 
squares:our data, full triangles: ref. J|) and icosahedron (open circles) at 
C(L)/L = 0.75 



• The continuum limit may very well be again higher than the prediction 
of Caracciolo et al (1.8810(3)), which most likely comes from a tran- 
sient, PT dominated regime; with present day's computers one cannot 
study reliably the larger L values needed to settle this issue. 

Our conclusion is that there is no evidence for the crossover phenomenon 
in the discrete nonabelian symmetry enhancement suggested by Caracciolo et 
al. While it follows from the accepted PT calculations, their scenario seems 
bizarre from the point of view of the discrete spin model. Indeed the latter 
may simulate a continuous symmetry only if the spins are correlated over a 
sufficiently large portion of the lattice. Fluctuations alone are not sufficient 
to enhance the symmetry, as can be seen by the high temperature expansion, 
which clearly is different for the icosahedron and the 0(3) models. Symmetry 
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enhancement requires the collective effect of many discrete spins over the 
typical distance of a correlation length; therefore one should expect that 
the discrete spins are capable of approximating the continuous ones better 
and better with increasing correlation length. Thus it would be bizarre if 
as the correlation length increases, beyond some value, correlation functions 
at distances measured in units of the correlation length would become less 
symmetric. However, being based only on some numerics, the present work 
cannot rule out the Caracciolo et al scenario. 

What the present work does show though is that the scaling curve pre- 
dicted by Caracciolo et al in 1995 [[J is not correct and that the lattice 
artefacts are larger than they claimed. If one accepts the universality of the 
0(3) and icosahedron spin model, then it appears that one could get a much 
better approximation of the continuum scaling curve by studying the icosa- 
hedron spin model. For instance the scaling curve produced by Caracciolo et 
al never reached 2, presumably because the best it could do is reach the 
value predicted by perturbation theory. Now in the icosahedron model one 
can work at f3 crt . The results, shown in Tab. 7 show that £(2L)/£(L) does 
actually reach 2 (definition of (3 crt ) and that this happens at x ~ 1 (non- 
trivial information). If the scaling curve were known, one could repeat the 
Kim- Caracciolo et al procedure and produce some thermodynamic values for 
£(/3). It should be emphasized though that with that procedure, the error 
compounds and one can lose control very fast. For instance if the scaling 
curve predicted by Caracciolo et al is changed by adding to their fit for- 
mula for £(2L)/£(L) the term 0.02a; and the starting values were L = 100 
£ = 75 then the predicted thermodynamic value changes from 2000 to 2670. 
Please notice that the modification we introduced vanishes for x — > and it 
amounts to an increase of less than 1% for any x, yet the change in the pre- 
dicted value is about 33%. Since in their 1995 paper Caracciolo et al stated 
themselves that they did observe scaling violations of about 1.5% (and the 
data shown in Fig.l suggest that that was a gross underestimate), we believe 
that their claim of having verified AF at correlation lenghts 10 5 at the 4% 
level is unjustified. 



Note Added in Proof: 

After the completion of this paper, S. Caracciolo drew our attention to 
some recent work done by himself with A. Montanari and A. Pelissetto con- 
cerning the issue of nonabelian discrete symmetry enhancement 0, |14| . In 



this work Caracciolo et al claim to have shown that in fact the continuum 
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limit of the model perturbed by a term enjoying only icosadral symmetry is 
different from that of the 0(3) model; but in contrast to their earlier estimate 
that the difference should show up at correlation length of about 200 they 
now say that they have to go to correlation length of about 10 5 . 

Unfortunately the authors base their conclusion only on data taken on 
relatively small lattices (while their estimate of the thermodynamic correla- 
tion length is ~ 10 5 , their largest lattice for which they determine the step 
scaling function only is 150 x 150. As we already stated both in previous 
papers 0,0,0 and repeated throughout the present paper, such data can- 
not be considered as representing the true continuum limit since they place 
the system in a perturbative regime, known rigorously not to occur in the 
true continuum limit. It is also not clear if the perturbed system with the 
icosahedral symmetry is not in its magnetized phase at the values of the 
paramters chosen (the authors do not provide any infromation on this). 

The data reported in the present paper, at x — 0.25 and x = 0.5, do 
not suffer from this limitation and are consistent with discrete symmetry 
enhancement. In support of their claim, Caracciolo et al state that their data 
indicate that the difference between the model enjoying a discrete symmetry 
and 0(3) is increasing with the size of the lattice L. We have two comments 
to this observation: 

• Even though the continuum limit of the two models may be the same, 
there is no guarantee that this limit must be reached in a monotonic 
fashion, hence the observed increased descrepancy with increased L 
could be a transient phenomenon. 



• In some cases (Fig. 2 of [13|, Figs. 4 and 6 of ||14|| ) the data at their 
largest L value L = 150 seem to lie below the data at L = 90, contrary 
to their claim. 

We thank S. Caracciolo and A. Sokal for correspondence and discussions. 
A. P. wishes to thank the Werner Heisenberg Institut, where this study was 
conducted, for its hospitality. 
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Tab.l: 

£(L) and £(2L) /or the 0(3) mode/ a« £(£)/£ « .25 



L 


20 


40 


80 


160 


320 




1.332 


1.488 


1.6135 


1.728 


1.838 


£(£) 

e(2L) 


5.011(3) 
5.164(3) 


10.070(7) 
10.395(5) 


19.919(11) 
20.529(10) 


39.881(25) 
41.097(26) 


79.988(48) 
82.458(47) 



Tab.2: 

£(L) and £(2L) for the 0(3) model at £(L)/L w .5 



L 


20 


40 


80 


160 


320 


640 


(3 


1.595 


1.7143 


1.825 


1.935 


2.047 


2.16 




10.007(9) 
15.205(13) 


19.999(17) 
30.629(27) 


40.002(28) 
61.451(41) 


79.884(52) 
123.079(91) 


159.73(13) 
246.99(26) 


319.29(37) 
493.94(68) 



Tab.3: 

£(L) and f (2L) for the 0(3) model at £(L)/L w .75 



L 


20 


40 


640 


(3 


2.092 


2.18 


2.63 




15.010(16) 


29.607(30) 


470.07(86) 


i{2L) 


28.260(32) 


55.608(66) 


882.8(2.5) 



Tab.4: 



£(L) and £(2L) for the icosahedron model at £(L)/L ~ .25 



L 


20 


40 


80 


160 


320 


(3 


1.319 


1.457 


1.556 


1.6295 


1.6836 


Z(L) 
i{2L) 


5.016(3) 
5.164(3) 


10.043(6) 
10.348(4) 


19.960(14) 
20.583(11) 


39.955(15) 
41.237(13) 


80.278(57) 
82.829(41) 



Tab.5: 

£(L) and £(2L) for the icosahedron model at £(L)/L fa .5 



L 


20 


40 


80 


160 


320 


640 


(3 


1.5382 


1.618 


1.6767 


1.718 


1.7469 


1.7664 


i{2L) 


9.991(9) 
14.955(14) 


19.996(17) 
30.141(20) 


40.057(25) 
61.076(42) 


80.131(74) 
122.58(13) 


160.04(20) 
246.14(35) 


318.62(32) 
487.17(70) 



Tab.6: 

£(L) and £(2L) for the icosahedron model at £(L)/L fa .75 
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L 


20 


40 


80 


160 


320 


640 




1.733 


1.758 


1.774 


1.788 


1.794 


1.797 


£(2L) 


15.114(31) 
27.279(43) 


29.936(68) 
54.227(083) 


58.99(11) 
107.69(17) 


121.20(21) 
222.97(33) 


242.68(39) 
448.06(73) 


479.40(5) 
882.2(1.6) 



Tab. 7a: 

£(L) for the icosahedron model at (3 = 1.802 



L 


20 


40 


80 


160 


320 




19.946(49) 


39.50(11) 


78.11(21) 


155.77(38) 


305.20(51) 



Tab. 7b: 

£(L) for the icosahedron model at j3 = 1.803 



L 


320 


640 




318.23(76) 


642.2(1.4) 
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